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Abstract: We investigate the indirect effects of heavy vector bileptons being exchanged 
in polarized M0ller scattering, at the next generation of linear colliders. Considering both 
longitudinal and transverse beam polarization, and accounting for initial state radiation, 
beamstrahlung and beam energy spread, we discuss how angular distributions and asymme- 
tries can be used to detect clear signals of virtual bileptons, and the possibility of distin- 
guishing theoretical models that incorporate these exotic particles. We then estimate 95% 
C. L. bounds on the mass of these vector bileptons and on their couplings to electrons. 
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1 Introduction 



Several schemes have been put forward to address important questions left unanswered in 
the standard model. There are models in which some of the standard model fundamental 
particles, such as quarks and leptons, are thought to be in fact composite. Another approach 
is to extend the symmetry by enlarging the local gauge group. In an effort to incorporate 
gravity, more exotic proposals have been made, such as those which depend on the existence 
of extra spatial dimensions. Since the standard model has passed severe experimental tests 
at the e"'"e~ and hadron colliders, showing its validity up to the electroweak scale, one 
expects that it could be seen as an effective theory, valid up to some large mass scale A. 
Most theoretical arguments point to A ~ 1 TeV. Extended electroweak models predict the 
existence of new particles and interactions. An interesting class of exotic particles known 
as bileptons [Tj is present in several of these models, such as left-right symmetric models, 
technicolor and theories of grand unification. The bileptons in which we are interested 
are vector bosons which couple to standard model leptons and carry two units of lepton 
number. In particular, heavy gauge bileptons may be found in models where the standard 
SU{2)l X U{1) group is embedded in a larger gauge group. This is the case of some of the 
331 gauge models [2j. In this paper we propose to look for indirect signals of doubly-charged 
vector bileptons in M0ller scattering, with polarized beams, and at the new linear collider 
energies. The next generation of e^e~ colliders will be designed to operate at center-of-mass 
energies ranging from the Z° mass to a few TeV, with very high luminosities. Electron 
and positron beam polarizations are expected to be available at these facilities, providing 
experimentalists with powerful tools to carry out precision tests of the standard model 
and to explore new physical phenomena. Currently, the best lower bound on the vector 
bilepton mass is My > 850 GeV , a result which was established from muonium-antimuonium 
conversion [5] . Another useful lower bound My > 740 GeV has been derived [3] from current 
experimental limits on fermion pair production at LEP and lepton-flavor violating charged 
lepton decays. While less stringent, this limit does not depend on the assumption that the 
bilepton coupling is flavor- diagonal. 

Studies of bilepton effects in M0ller scattering have been made before [5], mostly in the 
context of SU(15) grand unified theories. In this paper we work in the framework of the 331 
models, concentrating on the minimal version, and extend the previous studies by taking 
into account (i) important beam effects such as initial-state radiation, beamstrahlung and 
beam energy spread; (ii) longitudinal and transverse polarization of the colliding beams, 
which are expected to be available in the next generation of linear colliders; (iii) Gaussian 
smearing of the four-momenta of the final-state leptons, simulating the uncertainties in the 
energy measurements in the electromagnetic calorimeters. All these items are significant 
for a realistic comparison with experimental data. We analyze several distributions and 
asymmetries in M0ller scattering at the Next Linear Collider (NLC) energies, with the 
purpose of searching for indirect signals of gauge bileptons. The predictions of an SU{lb) 
GUT model [6] for these observables are also shown for comparison. From the angular 
distribution of final-state electrons we establish bounds on the couplings and mass of these 
bileptons at the 95% confidence level. If evidence of the existence of these vector bileptons 
is found in M0ller scattering, one can also verify the validity of relations between the masses 
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of the vector bileptons and new neutral gauge bosons, which are connected to the Higgs 
structure of the 331 models. In section II we give a brief review of the 331 models [7], 
outlining only those features that are most relevant to our analysis. Section III describes 
in detail the numerical simulation of the M0ller events and the corresponding analysis of 
the results. 95% confidence level limits are established in section IV. Our conclusions are 
summarized in section V. 

2 Review of the 331 models 

The 331 models are based on the SU{3)c ® SU{3)l ® U{l)x gauge group and predict new 
physics at the TeV energy scale. They arrange the ordinary leptons in SU{3) l antitriplets, 
two generations of quarks in triplets and a third generation in an antitriplet. Anomaly 
cancelation, essential to a gauge theory, takes place not within each generation but when all 
three families of quarks and leptons are taken together. The fact that one quark generation 
transforms differently from the other two is essential in these models, since the number of 
triplets must equal the number of antitriplets to ensure that the models remain anomaly-free. 
This implies that the number of generations is divisible by the number of colors. Besides 
standard gauge bosons, the models predict a new neutral gauge boson Z' and four vector 
bileptons and Y^^. These gauge bosons are expected to be considerably heavier than 
the standard gauge bosons. In addition to the ordinary quarks, the class of models we are 
considering here contains three new heavy quarks, with exotic electric charges 5/3 and —4/3. 
However it is possible to build 331 models in which these exotic quarks are not present [8], 
and as a consequence, neither are doubly- charged gauge bileptons. Following the notation of 
ref . [S] , the coupling of the doubly-charged vector bileptons F to electrons and positrons 
is given by the interaction Lagrangian 

Unt = -^Y^^e^Cn^e - ^Y;-eri,Ce^ ■ (1) 

For the minimal 331 model, symmetry breaking can be accomplished by three complex 
SU{3)l triplets and a complex sextet, allowing these nonstandard gauge bosons to acquire 
plausible masses. In this case, the mass Mz' of the neutral vector boson Z' and the mass 
My of the doubly-charged bilepton are related by 

My ^ ^/3il-Asm'ew) 
Mz> 2cos^vK 

Alternative Higgs structures are possible in different 331 models, but then the relation above 
no longer holds. 

3 Simulation and analysis 

One substantial source of beam energy degradation is the initial-state radiation (ISR). This 
is a QED effect, and corresponds to the emission of photons by the incoming electrons and 
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positrons. For a considerable fraction of events, this bremsstrahlung emission lowers the 
effective center-of-mass energy available for the hard scattering process. To account for ISR, 
we used the structure function approach discussed in [9]. In order to achieve high luminosi- 
ties at the new linear colliders, the beams must have very small transverse dimensions. The 
particles in a colliding bunch suffer considerable transverse acceleration due to the collective 
electromagnetic fields produced by the particles in the opposite colliding bunches, which 
gives rise to the emission of synchrotron radiation, the so-called beamstrahlung. The effec- 
tive energy available for the reaction is then smaller than the nominal value. The average 
energy loss of a positron or electron by beamstrahlung depends on the design parameters 
of the accelerator. For some designs of a 1 TeV NLC, for instance, the colliding beams 
may lose about 13% of the nominal energy from beamstrahlung emissions. This loss may 
reach 31% at a 3TeV CLIC |lOj. To obtain the beamstrahlung spectrum, we followed the 
approach of ref. [IT], which is based on the Yokoya-Chen approximate evolution equation 
for beamstrahlung [12]. In this paper, the calculations of the beamstrahlung energy spectra 
were carried out starting from the energy-dependent sets of NLC design parameters [13]. 
The spectra corresponding to ISR and beamstrahlung emissions were convoluted, and the 
resulting distribution was used to compute all the required differential cross sections. We 
also considered a possible beam energy spread, which was taken to be Gaussian distributed, 
with a width of 1% of the nominal beam energy. 

In addition to the standard Feynman diagrams for M0ller scattering, in 331 models the 
process also proceeds via an s-channel exchange of a doubly-charged bilepton, and a t-channel 
exchange of a Z' boson as well. By ignoring beam effects, neglecting the Z' exchanges and 
considering only unpolarized beams, we verified that our numerical calculations agree with 
the trace calculation of ref. [5]. Likewise, we cross-checked our calculations with those of 
ref. [13], which were carried out in the framework of the standard model, but with arbitrary 
beam polarization. 

The differential cross sections were calculated with Monte Carlo techniques, with the sim- 
ulated events selected according to the following set of cuts : (i) the final-state electrons 
and positrons were required to be produced within the angular range |cos^j| < 0.95, where 
6i stands for the polar angle of the final-state lepton with respect to the direction of the 
incoming electron beam; (ii) all events in which the acoUinearity angle ( of the final-state 
e"*" — e~ three-momenta did not pass the cut ( < 10° were rejected; (iii) the ratio of the ef- 
fective center-of-mass energy to the nominal center-of-mass energy for any acceptable event 
was required to be greater than 0.9 . Our simulations were made in the context of the NLC 
program [TH], considering a maximum center-of-mass energy of 1 TeV, and several years of 
operation, so as to accumulate an integrated luminosity of 500 fb~^. In order to simulate 
the finite resolution of the NLC electromagnetic calorimeters, we Gaussian-smeared the four- 
momenta of the produced electrons and positrons [16] . The energies of the final-state leptons 
were distributed as a Gaussian with half-width AE of the form AE/E = 10%/\/E © 1%, 
whereas the directions of the lepton three-momenta were smeared in a cone around the cor- 
responding original directions, whose half-angle is a Gaussian with half-width equal to 10 
mrad. 

Beam polarization will play a useful role at the new linear colliders [17] . By using polarized 
lepton beams, one can effectively reduce backgrounds and increase the sensitivity of spin- 
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dependent observables to potential new physics. In our calculations we worked with the 
electron beam projection operators in the extreme relativistic regime, 

V{pa,b) = lim \{^a,b + m){l + 75 

^(1 + Ss) A,b + ^PT%5(cOS0a,fe f^l + sin(j)a,b ^2) ^a,b , 

where n'^f^ represent the spin vectors and P^'^(P^'^) stand for the longitudinal (transverse) 
polarizations of the incoming electron beams, whose four-momenta for a nominal center-of- 
mass energy y/s are given by = (^, 0, 0, ^) and p'^ = (^, 0, 0, — ^). For the numerical 
calculations dealing with transverse polarization of the electron beams, the azimuthal angles 
of the transverse polarization vectors were taken to be (pa = 4'b = 0, and the purely spatial 
vectors < = (0, 1, 0, 0) and = (0, 0, 1, 0). 

In order to discuss how to probe the presence of gauge bileptons in M0ller scattering, we 
first considered longitudinally polarized electron beams. For the longitudinal polarizations 
of the beams we considered P£ = —Pi = 0.9, with an uncertainty given by APl/Pl = 0.5%. 
Figure 1 shows the dependence of the total cross section on the center-of-mass energy, for 
an input bilepton mass My = 1.2TeV. Only for higher energy values is the total cross 
section significantly altered by the exchange of a vector bilepton. The angular distribution 
da/d{cos6) of the final-state electrons is, however, more sensitive to the presence of such 
a particle. This is displayed in Fig. 2, where the angular distribution is plotted both for 
^/s = 500GeV and ^/s = 1 TeV, with the corresponding curves for the standard model 
shown for comparison. The M0ller scattering angular distribution in the minimal 331 model 
differs from that of the standard model for most of the angular range, the more so for a 
center-of-mass energy ^/s = 1 TeV. These deviations from the standard model predictions 
can be used to establish bounds on the mass of the gauge bilepton, and its coupling to 
electrons. The sjnnmetric shape of the angular distribution suggests that the integrated 
forward-backward asymmetry Aps should be small. This is indeed the shown in 

Fig. 3, where Aps is plotted for several input values of My, at an energy a/s = ITeV, and 
the one-standard- deviation error bars represent only the statistical errors. Next we analyze 
the discovery potential of spin asymmetries. 

Starting from the polarization-dependent angular distributions, one can compute the follow- 
ing asymmetries: 

4/ n^ dO\^-\i p\,-\i -r uuy^-\i \i - uuy^\i p\) - uu^\i p\, \i 

/ill cosy I = — — ; — — , „, ,. ; — ; — — — , „, ,. ; — ,, „ , , „, ,. ; — ,, „ , — , „, ,. lol 

\ ' j_/ I nal I DoW I J_/ I r>a\ I r>b\\ < j — Dal I nbw i nal I nbw ^ ' 

(4) 

= d.(-|P£|.|Pl|)+MO,0) 

In the limit |P£| = |P^| = 1, A^ (cos9) and A2{cos9) reduce to the familiar parity violating 
M0ller asymmetries [TS] 
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respectively. 

The behavior of spin asymmetry Ai as a function of cos6 is depicted in Fig. 4, for a bilepton 
mass My = 1.2TeV. At \/i = SOOGeV^ the deviation from the standard model values 
is small. This asymmetry becomes more sensitive to the presence of a vector bilepton at 
an energy of 1 TeV, and its angular dependence distinguishes the two nonstandard models . 
Statistical errors for the spin asymmetries discussed in this section are rather small, and most 
of the effects of the systematic errors are expected to cancel out in these asymmetries. Fig. 5 
shows that a similar pattern holds for A2, which is represented as a function of cos6 for both 
500 GeV and 1 TeV. At this latter energy, however, A2 does not lead to a clear distinction 
between the minimal 331 model and the SU (15) GUT model. Asymmetry A^ concerns the 
difference between a polarized angular distribution and its unpolarized counterpart, and 
is displayed in Fig. 6. Unlike Ai and A2, which take values of the order of 5%, A3 may 
become quite large, reaching a maximum magnitude of about 63% for the minimal 331 model 
at y/s = 1 TeV. It can be useful to discriminate a model with vector bileptons from the 
standard model, even at a 500 GeV linear collider. At a center-of-mass energy of 1 TeV, the 
difference between the prediction of the 331 minimal model for A^ and the standard values 
is even more striking. 

We assume that in the next generation of linear colliders, transverse polarization of electron 
beams will be available as an extra tool to search for the new physics. This could be achieved 
by means of spin rotators, which convert longitudinal into transverse polarization. It is not 
immediately clear whether the use of transversely polarized electron beams may add any 
important information, which could not be extracted from longitudinally polarized M0ller 
scattering. In fact, standard M0ller scattering is known not to be particularly sensitive to 
transverse beam polarization [T3], but the presence of an s-channel vector bilepton might 
modify this picture. Effects of transverse beam polarization in the nonstandard M0ller 
scattering under discussion only materialize if both beams are transversely polarized. In 
order to search for possible advantages of transverse beam polarization for the problem at 
hand, we examined in detail the behavior of the following differential azimuthal asymmetry: 



J (+)"''+' d{cose)d4, J(- 
/(+) '^'^ d[cose)d<j> + /(-) d(cos9)dcj) 



where the subscript +(— ) indicates that the integration over the azimuthal angle (j) is to 
be carried out over the region of phase space where cos2(j) is positive (negative). As far as 
Ait{cos6) is concerned, and considering an integrated luminosity of 500 we found that 
it would be difficult to separate the signal from the standard model background at center- 
of-mass energies around 500GeV^, even if My is only moderately large. At ^/s = 1 TeV, 
however, it is feasible to detect the effects of a vector bilepton, as long as its mass is not 
much larger than the center-of-mass energy. This is illustrated in Fig. 7, where Ait{cos6) is 
plotted both for an input mass My = SOOGe^ and and for My = 1.2TeV, along with the 
standard model expectation. We also investigated the integrated version of the asymmetry 
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above, 



^2T — - — — -inz — — — ttt, — 12Z — ^ 



/(+) d{cose)d<p + d{cos9)d^ 

where the integrations are consistent with the cuts specified in section III. The mass de- 
pendence of A2T for ^/s — 500GeV and ^/s — 1 TeV is presented in Fig. 8. A2T is found 
to be sizable over a fairly wide range around the bilepton resonance. We checked that the 
shape of the curve representing the mass dependence of this asymmetry becomes wider as 
g^i increases, while the position of the corresponding minimum remains essentially the same, 
for a fixed center-of-mass energy. 



4 Bounds 

A test was apphed to estimate discovery hmits for a vector bilepton in M0ller scatter- 
ing. Considering only longitudinal beam polarization, we compared the angular distribution 

da / d{cos9) of the final-state electrons, modified by the presence of a vector bilepton, with 
the corresponding standard model distribution. Assuming that the experimental data will be 
well described by the standard model predictions, we defined a two-parameter estimator 

x\9.uMy)=Y,i^\^) (10) 
i=\ « 

where N^^ is the number of standard model events detected in the i*'* bin, Ni is the 
number of events in the i*'* bin as predicted by the model with bileptons, and ANf^ — 

JiVW^y + (A'f*^e)2 the corresponding total error, which combines in quadrature the 
Poisson-distributed statistical error with the systematic error. For the latter we assumed 
a conservative value e = 5% for each measurement. The angular range \cos9\ < 0.95 was 
divided into Nh = 20 equal-width bins. The coupling g^i of a vector bilepton to electrons and 
the bilepton mass My were varied as free parameters to determine the distribution. The 
95% confidence level bound corresponds to an increase of the by 5.99 with respect to the 
minimum Xmin the distibution. Fig. 9 presents the resulting 95% C. L. contour plots on 
the {g^i, My) plane for the nominal center-of-mass energies ^/s = 500GeV and y/s = 1 TeV. 
The unpolarized case is shown in Fig. 10. Since our differential cross sections contain only 
even powers of g^i, it suffices to use positive values of the coupling in Figs. 9-10. We also 
calculated the corresponding 95% C. L. limits on the bilepton mass at these NLC energies, 
considering only the minimal 331 model, in which a Z' exchange has to be taken into account. 
The results are displayed in Table I. 



5 Conclusions 

Future linear colliders will provide an opportunity to look for new particles and their in- 
teractions. In this paper we discussed how to search for effects of a virtual vector bilepton 



7 



Table 1: 95% C. L. limits on the bilepton mass in the minimal 331 model, at NLC energies 
Polarization ^^=500 GeV ^/s=l TeV 
unpolarized 1230 GeV 1815 GeV 
polarized 2529 GeV 4574 GeV 

in polarized M0ller scattering. Starting from a realistic simulation of this process, we ana- 
lyzed several polarization-dependent observables that might provide strong evidence of the 
existence of vector bileptons, and allow to discriminate models or classes of models which 
predict these particles, should any deviation from the standard expectations be detected. 
We demonstrated that these polarization-dependent observables are more sensitive to vector 
bileptons than their unpolarized counterparts, in a large region of the allowed parameter 
space. This sensitivity was found to be stronger at a center-of-mass energy of ITeV than 
at 5006*6^^. The bounds on the masses and couplings of the bileptons, derived from a 
estimator, indicate that it should be possible to probe mass scales of up to several TeV for a 
signal of a vector bilepton. Longitudinal polarization of the electron beams in M0ller scatter- 
ing has proved useful to improve these bounds on the masses and couplings. Although most 
of our calculations were carried out in the context of the minimal 331 model and an SU{15) 
GUT model, we believe that our overall conclusions could be extended to other models with 
bileptons. 
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Figure 2: Angular distribution of the final-state electrons; Curves (A) and (D) show the 
angular spectra predicted by the standard model at \/i = SOOGeV" and = 1 TeV respec- 
tively, while (B) and (C) display the corresponding angular distributions in the minimal 331 
model. Curve (E) shows the prediction for the SU(15) model at ^/s — ITeV. 
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Figure 3: Forward-backward asymmetry Aps for several input masses My, at y/s — ITeV, 
according to the minimal 331 model (B). The upper curve (A) shows the expected value for 
the standard model 
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Figure 4: Polar angle dependence of spin asymmetry Ai{cos6)] Curve (A) shows the pre- 
diction for the SU(15) model at i/i = ITeV. Curves (D) and (B) show the standard 
model predictions for Ai{cos9) at ^/s = bOOGeV and ^/s — ITeV respectively, while (C) 
and (E) represent the corresponding expectations for the minimal 331 model, for a mass 
My = 1.2TeV. 



13 



T 




Figure 5: Polar angle dependence of spin asymmetry ^2(005^^); Curve (A) shows the pre- 
diction for the SU(15) model at i/i = ITeV. Curves (B) and (D) show the standard 
model predictions for ^42(005^) at ^/s = bOOGeV and ^/s — ITeV respectively, while (C) 
and (E) represent the corresponding expectations for the minimal 331 model, for a mass 
My = 1.2TeV. 
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Figure 6: Polar angle dependence of spin asymmetry ^3(005^^); Curves (B) and (C) show 
the standard model predictions for 743(0056*) at \/i = 500GeV and = 1 TeV respectively, 
while (D) and (A) represent the corresponding expectations for the minimal 331 model, for a 
mass My = 1.2TeV. Curve (E) shows the prediction for the SU(15) model at = ITeV. 
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Figure 7: Transverse polarization asymmetry Ait{cos9) as a function of cos6, at = 1 TeV. 
The curve with error bars (D) corresponds to the minimal 331 model with My = SOOGeV, 
whereas curve (C) corresponds to the SU{15) GUT model prediction for the same mass. 
Curve (B) represents the behavior of the asymmetry for a 1.2 TeV bilepton in the minimal 
331 model, and the histogram (A) corresponds to the standard model expectation. The error 
bars for the upper curves are similar to those of the lower curve on a bin-to-bin basis. 
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Figure 8: Transverse polarization asymmetry A2T as a function of My. The solid histogram 
(A) refers to the minimal 331 model prediction at an energy ^/s = bOOGeV and the long- 
dashed histogram (B) to the SU{lb) model at ITeV. The resulting asymmetry A2T for the 
minimal 331 model at = 1 TeV is is represented by the dashed histogram (C). 
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